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We investigate frustrated antiferromagnetic Heisenberg 
quantum spin chains at T = for S = | and S — 2 us- 
ing the DMRG method. We localize disorder and Lifshitz 
points, confirming that quantum disorder points can be seen 
as quantum remnants of classical phase transitions. Both in 
the S = | and the S — 2 chain, we observe the disappearance 
of effectively free S = i and S = 1 end spins respectively. 
The frustrated spin chain is therefore a suitable system for 
clearly showing the existence of free end spins S = [S/2] also 
in half-integer antiferromagnetic spin chains with S > \. We 
suggest that the first order transition found for S — 1 in our 
previous workcl is present in all frustrated spin chains with 
S > i, characterized by the disappearance of effectively free 



end spins with S' = [S/2]. 
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I. INTRODUCTION 

Because of its particularly simple structure the 
isotropic antiferromagnetic (AFM) Heisenberg quantum 
spin chain with frustration through an AFM next- 
nearest-neighbor (NNN) interaction at T = 0, 



(«>0), (1) 



can be considered a textbook example for a system with 
competing interactions. In recent years, this model has 
been studied extensivelyn because of its theoretical in- 
terest and experimental realizations. However, so far 
mainly results for spin lengths S = h and S = 1 are 
availabletrQ. Although these cases are the most rele- 
vant ones for experiment, the underlying theory cannot 
be considered complete; e.g., the generalization to large 
S and the emergence of the classical limit are unclear. 
In this paper, we will first outline the available results, 
and take a step towards large-5 behavior by numerically 
investigating the S = § and S — 2 case. Our results, 
among them a clear demonstration of free end spins in 
the S = | chain, will allow us to extract the large-S" 
behavior and to propose a scenario for frustrated chains 
with arbitrary S accommodating available results. 



II. S = ~ AND S = 1 FRUSTRATED CHAINS 

Unfrustrated (a = 0) chains show jitrongly different 
behavior for half-integer and integer SB; the same holds 
if frustration is switched on. 

The frustrated spin-^ chain is characterized by a cr 
ical phase with AFM correlations for a < ac — 0.24111 
At ac, there is a continuous phase transition to a dimer- 
ized gapped phase for all a > ac- The dimerized phase 
is characterized by the non-vanishing dimerization order 
parameter 



D{a) = |(Si-S 



i+l 
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which reaches its maximum for a — 0.5781 and vanishes 
exponentially for large aa. The gap opens exponentiallyu, 
reaches a maximum at a « 0.6, and disappears in the 
large-cv limit, where the frustrated chain decomposes 
into two unfrustrated critical Heisenberg chainsQ. This 
scenarioJJlustrates the generalized Lieb-Schultz-Mattis 
theoremt3, which allows just these two types of phases 
for a half-integer spin chain with rotational and trans- 
lational invariance. While the phase transition does not 
affect the rotational invariance of the ground state, par- 
ity is broken. For S = |, the transition is suggested to 
be at ac ~ 0.33, and it can be reasonably concluded that 
it generalizes to all half-integer spinaLi. 

The frustrated spin-1 chainlil is gapped for all a. For 
a < ac — 0.745, one finds a phase which can be char- 
acterized by the S = 1 Affieck-Kennedy-Lieb-Tasaki 
(AKLT) madelij. This phase has non-vanishing nonlocal 
string orderE3 and effectively free S = | spins at the ends 
of open chains. These give rise to a low-lying excitation 
triplet, which is degenerate .arith the ground state in the 
thermodynamic (TD) limito. While this excitation is 
clearly a boundary effect, it is linked to bulk behavior, as 
shown by the existence of nonlocal string order. At ac, 
there is a first order transition characterized by discon- 
tinuous disappearance of the nonlocal string order and 
the edge excitation triplettl. The large-a phase without 
string order or free end spins can be understood by a 
next-nearest neighbor generalization of the AKLT-model 
which can be extended to all Sid. 

Both chains exhibit disorder points and Lifshitz points 
(i.e. points, where the correlation function in real space 
becomes incommensurate, or where the structure func- 
tion develops a two-peak structure, respectively; for de- 
tails of the concept, see Rcfs. [ljjl5| and Refs. cited therein. 
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Some authors merely call the tricritical point, where clas- 
sical disorder and Lifshitz lines meet, a Lifshitz point.). 
For 5 = i, the disorder point ajj — 0.5 is the well- 
known Majumdar-Ghosh point, the Lifshitz point is at 
at, = 0.5206O. The respective results, for S = 1 are 
a D = 0.284(1) and a L = 0.3725(25)1. We interpret 
these as proposed in Refs. p]Jl5|: In the classical frus- 
trated chain, there is a T = phase transition from com- 
mensurate AFM order to incommensurate spiral order 
for a c i = 0.25. At finite T, the classical chain is dis- 
ordered, but should exhibit disorder and Lifshitz lines, 
in line with the usual explanation of classical disorder 
linealj. The link between classical chains at T ^ and 
quantum chains at T = is established by the nonlinear 
cr-model relation T oc 1/5, which holds at least for small 
frustration. This scenario has three implications: First, 
disorder and Lifshitz points are nothing but the quantum 
remnants of the classical phase transition, second, they 
should thus be present for arbitrary 5 < oo; and third, 
in the classical limit, their ppsitions should converge to- 
wards each other and to a c ;EJ. These implications remain 
to be tested. Let it be remarked that as a side result, 
we conclude that the transitions in the 5 = \ and the 
5 = 1 chain are pure quantum phenomena unrelated to 
the classical transition. _ . 

We apply now the DMRG algorithmtll to frustrated 
5=| and 5 = 2 quantum spin chains, considering chains 
of length up to 300 sites, keeping up to 300 states. All 
chains considered have open boundary conditions. We 
map out the chain behavior with a. 

III. DISORDER POINTS AND LIFSHITZ POINTS 

Numerically, we can identify disorder point and Lif- 
shitz point in both cases investigated, using the criteria 
already employed in Ref . [l| For 5 = § , we find a disorder 
point an — 0.386(1) and a Lifshitz point a_L = 0.388(1); 
they are close, but separable. The correlation length 
has a minimum £(«£>) ~ 1-9. For 5 = 2, disorder 
point and Lifshitz point are found for aj, = 0.289(1) and 
(Xu = 0.325(5), with a minimum of correlation length 
£(a D ) « 1.8. n 

Haldane's now-well established conjectures predicts a 
different behavior of half- integer and integer unfrustrated 
Heisenberg spin chains. For reasons of continuity, one can 
expect that this leads to substantial differences in not too 
strongly frustrated half-integer and integer spin chains 
also. From this one can expect that the disorder and 
Lifshitz points of half-integer and integer spin chains take 
different tracks to the classical limit. This is supported 
by putting known results together (see Table 1). 

Half-integer and integer spin chains show clearly dif- 
ferent behavior. For half-integer spin chains we find, that 
the disorder and Lifshitz points, which were already rel- 
atively close together for 5 = -|, are just distinguishable 
in the case of 5 = | and much closer to a c ; than their 



5 = \ counterparts. At the same time, the disorder point 
correlation length increases. As the 5 = \ disorder point 
is the Majumdar-Ghash ground state, only neighboring 
spins are correlated!!^, so that the correlation length is 
not defined or zero. 

For the 5 = 2 chain we find that disorder and Lifshitz 
points move, as expected, closer to each other, but not as 
close as in the half-integer case, together with an increas- 
ing correlation length at the disorder point. However, the 
disorder point of the 5 = 2 chain is slightly further away 
from a c i = j than for the 5=1 chain. This small excep- 
tion is in our view due to the fact that our model of the 
behavior of disorder and Lifshitz points does not take 
into account other physical effects in the chains which 
may influence the exact locations, but only the differ- 
ence between half-integer and integer spin chains. Once 
disorder and Lifshitz points have moved close to the clas- 
sical limit, these additional effects may be strong enough 
to enforce deviations from the by then slow convergence. 
Summarizing, we take these results as strong support for 
our explanation of quantum disorder points and the im- 
plications listed above. 
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FIG. 1. Magnetization of the lowest-lying spin-1 (solid) 
and spin-2 excitations (dashed) for the frustrated S = § chain 
at q = 0.385 with L = 120. The spin-1 excitation is clearly a 
boundary excitation (Kennedy triplet), the spin-2 excitation 
combines a bulk excitation with the spin-1 edge excitation. 
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S = 2, a = 0.44 
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FIG. 2. Magnetization of the lowest-lying spin-1 excita- 
tion for the frustrated S — 2 chain at a = 0.44 for L = 200, 
below the point of the disappearance of edge excitations. 
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FIG. 3. Magnetization of the lowest-lying spin-1 excita- 
tion for the frustrated S = | chain at a = 0.60, which is a 
bulk excitation. 



IV. FREE END SPINS IN OPEN SPIN CHAINS 



An important observation in the frustrated 5 — 1 chain 
was the existence of free 5 = 5 end spins disappearing 
for sufficiently strong frustration. No such observation 
can be made in the 5 = | chain. We observe, for weak 
frustration, the existence of free s = | end spins in the 
5=1 chain, giving rise to a low-lying edge excitation 
triplet, well known as Kennedy triplet in 5 = 1 chains. 
Ngcj has conjectured the existence of free end spins in all 
unfrustrated open chains for 5 > 1. He derived a Berry 
phase contribution to the effective action of the nonlin- 
ear er-model, predicting free end spins S' = [5/2] for all 
5 > 1. In the critical phase, edge .and bulk excitations 
are hard to distinguish numericallycj; the presence of a 
gapped phase in the 5=1 frustrated spin chain makes 
the identification very easy due to the finite correlation 
length. As example, we show in Figure I the magnetiza- 
tion of the Sl otal — 1 state of the spin-1 end excitation 
triplet for a — 0.385 as well as the magnetization of the 
lowest spin-2 excitation, which is the boundary excita- 
tion plus a bulk excitation. Similar observations can be 
made for 5 = 2 (Fig. ^), where one finds free 5 = 1 
spins, but similar observations had already beep-,clearly 
obtained for 5 = 2 at the unfrustrated pointc3~E3. 
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FIG. 4. Magnetization of the lowest-lying spin-1 excita- 
tion for the frustrated S — 2 chain at a = 0.52 for L — 200 
(~ 10£), above the point of the disappearance of edge excita- 
tions. 

Furthermore, we observe that the free end spins disap- 
pear, as in the 5=1 case, both for 5 = | and 5 = 2. 
The lowest excitations shown in Figs. || and || are clear 
bulk excitations. 

As explanation, for this result, we invoke the model 
proposed by Ngt2l. His derivation inherently assumes lo- 
cal Neel order: in his Berry phase contribution, local 
Neel order leads to alternating signs in a sum over local 
phase contributions, which can be interpreted as an in- 
tegral of the first derivative of these contributions, such 
that merely the boundary terms survive. For continuity 
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reasons, we may also expect this assumption to hold for 
weak frustration. For sufficiently strong frustration, this 
assumption fails, as we will find short-ranged spiral cor- 
relation. The free end spins should then disappear for 
sufficiently strong frustration. 

It might also be quite instructive to consider the va- 
lence bond toy model of the dimerized ground states of 
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the S = § and S 



| frustrated chains as shown in Fig- 



ure |5|. This toy model cannot be constructed for S = h . 
It is basically a conventional AKLT ground state, where 
the odd number of singlet bonds emerging from each site 
is split such between right and left neighbor that the dif- 
ference between bond numbers is minimal, i.e. 1, arguing 
that the underlying interaction is non-dimerizcd. This 
gives alternating "strong" and "weak" bonds. In an open 
chain of even length, there will be "strong bonds" at each 
end for minimal total energy, which leaves effectively free 
spins S' = [S/2]. For this model, an argumentation as 
for integer spins can be repeated to predict a first order 
transition for all half- integer S > | . Note that free £nd 
spins also exist if resonating bonds restore criticalitj£3. 
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FIG. 5. Valence bond toy model of the dimerized ground 
states of the S = § and S = | frustrated chains. Open cirles 
are sites; full circles spin-i fully symmetrized on each site. 
Straight lines are singlet bonds. Note the presence of free 
S — \ and S = 1 spins at each end. 

As in the S ^ 1 case the string order parameter is not 
characterizing the disappearance of the free end spins, 
the precise localization of the point of disappearance is 
difficult: in a small region, mixed bulk-boundary excita- 
tions will be energetically most favorable in finite length 
systems. Considering the lowest excitation gap, we were 
able to narrow down the points to a — 0.48(2) for S = | 
and a = 0.46(1) for 5 = 2. 

It is reasonable to assume that the observed scenario 
holds for all larger values of S also. 



V. CRITICAL AND MASSIVE PHASES 
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We want to show with our DMRG results, that the gap 
of the 5 = § chain opens according to (||) and use (|^) to 
determine a c | . Before we present our results, let us first 
point out, that because DMRG is a variational method 
and hence overestimates energies and energy gaps, fit- 
ting DMRG data with (|J) will return a a c 3 too small. 
Furthermore we are restricted to a region a > a c j , for 
which the gap is significant bigger than the numerical 
error. To test our procedure, we calculate the spin gap 
of the frustrated S = h chain close to its transition and 
used the known value of a c i — 0.24110. We find good 
agreement. 

To determine the gap of the S = § chain, we calculate 
the energies of the ground state and the first excitation 
for chains of up to L = 300 spins with up to M = 300 
block states kept in the range a = 0.35 ... 0.385. We 
extrapolate M — > oo and L — > oo to get the gap in the 
thermodynamic limit, a — 0.35 is the smallest value of a 
for which we can see a significant gap. The gap changes 
its behavior at approximately a = 0.375, so we fit our 
gap data in the range 0.35 < a < 0.37 with (||) and 

determine a r 3 0.29, a value that we consider as lower 

2 

bound of the true value. This is consistent with an earlier 
approximate resulttil a c 3 w 0.33. 




As for S = i, S = I has a Kosterlitz-Thouless phase 
transition from an antifcrromagnetic critical phase to a 
dimerized one, whereas the S — 2 chain is massive for all 
values of frustration. 

Therefore it is expected, that the spin gap in the | 
chain shows close tq-,the transition the same behavior as 
the i chain, namelyta 



FIG. 6. Onset of the gap in the S = | quantum spin chain 
at the Kosterlitz-Thouless transition. 

In Figure |^ we show the opening of the gap. Our 
DMRG results are represented by dots, the solid line 
shows the fit. It clearly can be seen that for a > 0.37 the 
gap changes its behavior. The smallest gap value we have 
is overestimated due to systematic DMRG errors. For 
large frustration, the gap decreases again monotonically. 
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The gapped phase is characterized by finite dimerization 
(Figure [?]), which reaches its maximum at a = 0.415(5), 
for smaller frustration than in the 5 — | case, and then 
also drops monotonically. While bigger in absolute terms 
than for S = h, relative to the ground state energy it is 
weaker, as it should be: in the classical limit, zero dimer- 
ization is expected, if the 5 2 factor governing the ground 
state energy is scaled out. 
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FIG. 7. Dimerization for the frustrated S - 



| spin chain. 



For 5 = 2, we have determined the correlation length, 
which shows that the chain is always massive. For large 
frustration, we expect that £ converges to 2£ (a = 0); 
the saturation effect in Figure || is purely artificial. It 
is well known that the DMRG underestimates correla- 
tion lengths, the more so, if the DMRG truncation error 
is large. For large frustration, the chain starts to de- 
compose into two chains, making the DMRG description 
worse, dramatically increasing the truncation error. The 
saturation effect comes from a growing underestimation 
of a growing correlation length. 
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FIG. 8. Correlation lengths for the frustrated S = 2 spin 
chain. Values are systematically underestimated for large £ 
by the DMRG. The saturation for large a is artificial, a slow 
convergence to a « 100 for a — » oo being offset by stronger 
underestimation due to a strong increase in the truncation 
error. 

The remaining important question is whether we can 
associate a first order transition with the disappearance 
of the free end spins as in the 5 = 1 case. We may 
expect that this transition, which is a pure quantum ef- 
fect, is softened for larger 5. Moreover, we do not have 
a clear indicator like the string order parameter in the 
5=1 case. This, together with the increasing numer- 
ical difficulties for larger 5, made it impossible for us 
to definitely answer this question. As the spin correla- 
tion length is finite for all spin lengths considered at the 
point of end spin disappearance, at least a second-order 
transition can be excluded. 

At the disappearance of the end spins, we observe a 
drop of the truncation error for a fixed number of block 
states by more than an order of magnitude, as opposed 
to much smoother change of the truncation error with a 
for all other values of frustration. This might be an in- 
dicator that a first-order phase transition has happened. 
Furthermore, just below the transition, spin-spin corre- 
lations with large spin-spin distances change drastically 
in character from bulk behavior. As we calculate corre- 
lations for spins symmetric around the chain center, the 
spins involved are close to the chain ends, indicating a 
coexistence phenomenon or a mctastability problem of 
DMRG. After the disappearance of end spins, the effect 
disappears. All this is however highly indirect evidence 
only, but would be consistent with a first-order transi- 
tion. Also for systematic reasons, we suggest that there 
is such a transition for all 5 > 1, but further work will 
be necessary to attack this question. 
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VI. CONCLUSION 



In conclusion, a coherent picture of the behavior 
of AFM quantum spin chains with NNN frustration 
emerges. The wealth of phenomena observed is due 
to a complex interplay of different features, only some 
of which are sensitive to the difference between half- 
integer and integer spins. Pure quantum phenomena co- 
exist with traces of the classical limit. The phase tran- 
sition from AFM to spiral order in the classical limit 
leaves quantum remnants in the form of disorder and Lif- 
shitz points, which converge towards the classical phase 
transition point in the large- S limit. This "classical" 
phenomenon is common to all half-integer and integer 
spins. All transitions in the quantum frustrated chains 
are therefore pure quantum phenomena and unrelated to 
the classical limit. Haldane's conjecture shows up in the 
different fate of unfrustrated half-integer and integer spin 
chains, when the frustrating interaction is switched on: 
whereas integer spin chains remain gapped for all values 
of frustration, half-integer spin chains exhibit a dimeris- 
ing continuous phase transition from a critical antifer- 
romagnetic phase to a dimerized gapped phase, a nice 
illustration of the generalized Lieb-Schultz-Mattis theo- 
rem. The third prominent feature is the disappearance of 
effectively free 5" = [5/2] end spins which are present in 
all these chains with S > h for low frustration. This gen- 
eralizes the previously found result in the S = 1 chain. 
Free end spins are thus also present in half-integer open 
AFM spin chains. It was not possible to clearly show that 
there is also a first-order transition present for S > 1. 
While it would be systematic to expect this, it should be 
kept in mind that the 5=1 Heisenberg chain, among 
integer spin chains, might be rather speciala. 

Most numerical calculations were carried out on a 200 
MHz PcntiumPro running under Linux. 



Spin length 


Old 






S = 1 

s = l 

S = 2 


0.5 
0.284(1) 
0.386(1) 
0.289(1) 


0.52063(6)0 
0.3725(25) 
0.388(1) 
0.325(5) 


1.20(2) 
1.88(2) 
1.78(2) 



TABLE I. Location of disorder points an and Lifshitz 
points cil, and the disorder point correlation length £d for 
various spin lengths. 
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